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Abstract 

The cratering history of main belt asteroid (2867) Steins has been in- 
vestigated using OSIRIS imagery acquired during the Rosetta fryby that 
took place on the 5 th of September 2008. For this purpose, we applied cur- 
rent models describing the formation and evolution of main belt asteroids, 
that provide the rate and velocity distributions of impactors. These mod- 
els coupled with appropriate crater scaling laws, allow the cratering history 
to be estimated. Hence, we derive Steins' cratering retention age, namely 
the time lapsed since its formation or global surface reset. We also investi- 
gate the influence of various factors -like bulk structure and crater erasing- 
on the estimated age, which spans from a few hundred Myrs to more than 
1 Gyr, depending on the adopted scaling law and asteroid physical parame- 
ters. Moreover, a marked lack of craters smaller than about 0.6 km has been 
found and interpreted as a result of a peculiar evolution of Steins cratering 
record, possibly related either to the formation of the 2.1 km wide impact 
crater near the south pole or to YORP reshaping. 
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1. Introduction 



The European Space Agency's (ESA) Rosetta spacecraft passed by the 
main belt asteroid (2867) Steins with a relative velocity of 8.6 km/s on 5 
September 2008 at 18:38:20 UTC. The Rosetta-Steins distance at closest 
approach (CA) was 803 km. During the flyby the solar phase angle (sun- 
object-observer) decreased from the initial 38° to a minimum of 0.27° two 
minutes before CA and increased again to 51° at CA, to reach 141° when 
the observations were stopped. A total of 551 images were obtained by the 
scientific camera system OSIRIS, which consists of two image rs: the Wide 



Angle Camera (WAC) and the Narrow Angle Camera (NAC) (IKeller et al. 



20071 ). The best resolution at CA corresponded to a scale of 80 m/px at the 



asteroid surface. 

Steins has an orbital semi-major axis of about 2.36 AU, an eccentricity of 
0.15 and an inclination of 9.9°. It is therefore orbiting in a relatively quiet re- 
gion of the main belt, far from the main escape gateways, namely the secular 
Vq and mean motion 3:1 resonances. Its shape can be fi tted by an ellipsoid 



having axis of 6.67 x 5.81 x 4.47 km (IKeller et all 120101 ) 



Previous space missions have visited and acquired detailed data for a to- 
tal of 5 aster oids, namely three main belt asteroids (951 Ga spra, 243 Ida, 



253 Mathilde; IVeverka et all Il999al : iBelton et all 11992. 119941) and two near 



Earth objects (433 Eros, 25143 Itokawa; IVeverka et all Il999bl : ISaito et al. 
20061 ). Itokawa is the smallest of them, with dimensions of 0.45 x 0.29 x 



0.21 km. The other asteroids have average sizes ranging from about 12 km 
to 53 km. In this respect, Steins with its 5.3 km size lies between Itokawa and 
the large asteroids. It is therefore the smallest main belt asteroid ever imaged 
by a spacecraft (except for Ida's satellite Dactyl with a diameter of roughly a 
km). Moreover, Steins is a member of the relatively rare E-type class (com- 
posed by igneous materials), while other asteroids visited by spacecraft are 
members of the most common S- and C-type classes. Previous spacecraft 
observations opened a new field of investigation, namely the cratering of as- 
teroidal surfaces. A number of interesting processes were therefore studied 
with unprecedented detail, like the cratering on low gravity bodies, regolith 
formation, seismic shaking (e.g. IChapmanl . 120021 ). 

This paper analyzes some of the highest resolution OSIRIS images with the 
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aim to study the crater size distribution and derive the chronology of the 
impacts on the surface of the asteroid. This will also provide clues on the 
Steins bulk structure, evolution, and give new insights on the above men- 
tioned processes. 



2. Steins crater population and geological assessment 



A total of 42 crater-like features with dimensions ranging from about 0.2 
to 2.1 km have been identified on three WAC images obtained near closest 
approach and one NAC image obtained about 10 min before CA (see fig. [I]). 
However, the illumination conditions of the WAC images were considerably 
different from those of the NAC. Therefore in the present work, we restrict 
the analysis to the WAC images only. The effective area over which counts 
have been performed is 23.7 km 2 , or approximately 25% of the estimated 
total Steins' surf ace of 9 7.6 km 2 (for more details on crater identification and 
size estimate see iBessd . |2010[ ). 

A total of 29 crater-like features larger than 3 pixels, ranging from about 
0.24 km to 2.1 km, have been identified and used for the present work. The 
largest crater (named Ruby hereinafter) has a diameter of 2.1 km. The true 
nature of some features is uncertain. In particular, this is the case for a 
feature consisting of 8 pits aligned along a straight chain crossing a large 
degraded crater and extending almost from the south to the north pole of 
the asteroid. Owing to lack of sufficient resolution, it is not clear whether or 
not these aligned pits are due to impacts. The probability for such a chain of 
impacts to occur on a low gravity body like Steins is indeed extremely low. 
Moreover, a NAC image (see fig. [1] bottom panel), obtained from a different 
aspect, shows the presence of a large depression and an array of rimless cir - 



cular pits, continuing the pit-chain imaged by the WAC ( Keller et all |2010| ) 



Therefore, the series of aligned depressions are most likely not caused by 
individual impacts. An alternative explanation is that the impact forming 
the crater Ruby triggered the nucleation of a long fracture or fault, whose 
expression at the surface is the pit-chain. S imilar linear featur es were also 
observed on other small bodies, like Gaspra (IBelton et all Il992h . Other un- 
certain crate r-like features are located close to the rim of the Ruby crater 
flBessel . l2O10h . 

The cumulative distributions of all detected crater-like features and 18 bona 
fide craters are reported in fig. |2j The latter has been constructed reject- 
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ing all the uncertain craters described above. Note that the small bona fide 
craters (D < 0.5 — 0.6 km) are strongly underrepresented with respect to the 
distribution of the crater-like features. 

The population of craters shows a wide range of depth-to-diameter ratios, 
varying from very shallow craters (~ 0.05) to deep cr aters ( y 0.2 5). The 
average crater depth-to-diameter ratio is 0.12 ± 0.05 flBessel . l2010h . Small 
craters have typically a lower depth-to-diameter ratio. These characteristics 
are consistent with crater degradation due to ejecta blanketing and/or dis- 
turbance of loose reRolit h on the surface triggered by impact se ismic shaking 



(|Richardson et al.l . l2005l ). in agreement with results for Itokawa ( jHirata et al. . 
2009h . Furthermore, Ruby has sharper rims and a higher depth-to-diameter 



ratio (> 0.14; the bottom of the crater is not visible in the images) than 
several impact craters. Concerning the size distribution of Steins' surface 
material, the modeling of the observed strong opposition effect and the over- 
all photometric properties suggest that Steins may be covered by a layer of 
regolith having a mean particle size of ~ 100 fim. In the light of these con- 
siderations, the observed large degree of crater degradation may be explained 
in terms of regolith blanketing, i.e. deepest and sharpest craters are younger 
than more degraded ones. 



3. The impactor flux 

The flux of impactors can be expressed in terms of a differential distri- 
bution, <f)(d,v), which represents the number of incoming bodies per unit of 
impactor size (d), impact velocity (v) and per unit time. Such a differential 
distribution, under the assumption that the impact velocity is independent 
on the impactor size3, can be written as: 



<f>{d,v) = h{d)f(v) (1) 

where h(d) is the impactor differential size distribution, and f(v) is the 

i mpact probabil i ty per unit i mpact veloc- 



distribution of impact velocity (i.e. 
ity) normalized to J f(v) dv = 1 ( 



sec 



Marchi et all . I2005L 120091 . for details). 



1 This assumption stems from the lack of systematic differences, at least for the available 
observational data, in the orbits of MBAs according to their sizes. 
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In principle, the impactor size distribution could be derived from observa- 
tions. At the small size range relevant for Steins (d < 1 km), however, there 
is too little observational information to enable such an approach. To over- 
come such lack of informatio n, we use th e average size distribution of the 
main belt model derived by iBottke et al.l (j20Q5bl Jal) . The estimated num- 
ber of impactors at Steins is then obtained by multiplying the average size 
distribution by the intrinsic probability of collision with Steins, Pi (see fig. 
|3] upper panel). The latter, evaluated by taking into account the observed 
orbital distribution of main belt asteroids, is Pi = 2.87 • 10~ 18 km _2 yr _1 , 
nearly equal to the main belt average intrinsic probability of collision. The 
impact velocity distribution f{y) for Steins has been evaluated considering 
the population of main belt asteroids that presently inte rsect its orbit (see 
fig. [S lower panel). Computations have been done using the lFarinella fc Davis 
(119921 ) algorithm. The resulting average impact velocity is 5.7 km/s. 



4. The scaling law 



The impactor flux is converted into a cumulative crater distribution (the 
so-called Model Production Function, MPF) using a scaling law (SL). The 
physics of the cratering processes for small asteroids is still poorly known. 
In t his work we use the most upd ated SL derived from experi mental anal- 



ysis (IHolsapple and Housenl . 120071 ) and hydro code simulations (INolan et al. 



19961 ). The SL for cratering on asteroids depends on several parameters, the 
most important being the internal structure and tensile strength (Y). Both 
are unknown for Steins. However, some constraints can be obtained from 
morphological studies. 



Let us first consider the IHolsapple and Housenl (120071 ) scaling laws (HSL, 
hereinafter). An interesting issue concerns the minimum impactor dimen- 
sion for catastrophic disruption (d c d) of Steins. Assuming that Steins is an 
unfractured silicate ro ck and using the relevant specific energy for disruption, 



Q* D = 1 - 2 • 10 7 erg/g flHolsappld - Eooih , we derive that an impact at an aver- 



age modulus velocity of 5.7 km/s with a body having size d c d = 0.20—0.25 km 
would be sufficient for catastrophic disruption. Moreover, for an unfractured 
rock with surface gravity g and density p, the transition f rom strength to 
gravi ty cratering occurs at a crater diameter of ~ 0.8Y/gp (lAsphaug et al. , 
19961 ) . which exceeds the size of the Ruby crater, except the cases of unrea- 
sonably low Y values for a silicate body (namely < 10 5 erg/g). Therefore, 
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using HSL the strength regime applies. Under these conditions, we obtain 
that a 2 km crater is produced by an impactor having d ~ d cc i- 
In this respect it is interesting to note that even though the visual appear- 
ance of Steins is dominated by the big Ruby crater, the ratio of largest crater 
diameter to average asteroid diameter of 0.38 is not particularly high. Aster- 
oids Id a (0.44), Mathi lde (0.62), and Vesta (0.87) reach considerably larger 
values (lAsphaugi . 120081 ) . However, if compared to the specific energy required 
to disrupt the body, the big crater of Steins stands out (Fig. @|. Therefore 
the existence of the Ruby crater may be an indirect evidence that Steins was 
not a solid rock at the time of Ruby formation. Even if a particular tuning 
of the parameters may leave open such a possibility, it is more likely that 
Steins was a rubble pile or a collection of cohesive rubble of rocks. In the first 
case, for a pure cohesiveless rubble pile the gravity scale would apply. In the 
latter cohesive rock scaling law may be more suitable. Concerning 

the present state of Steins, it is li kely a rubble pile independently of its state 



prior to the formation of Ruby (jJutzi et al.l . I201CH ) . This conclusion is also 



in agreement with the two large fracture-like features seen on its surface (see 
previous sections). A preliminary study of the for mation of Ruby cr ater has 



been recently performed via numerical modeling (jJutzi et al.l . |2010| ). It has 



been found that a rubble pile body with some micro-porosity would have sur- 
vived the formation of the Ruby crater, although the non-porous monolithic 
body hypothesis cannot be ruled out. 

A further indication in favor of the rubble pile (both cohesiveless or with 
some low cohesion) nature of S t eins m ay come from its shape, possibly due 
to YORP spin-up flKeller et all . l2010f ). 

In conclusion, from previous reasoning, we limit our investigations to HSL for 
cohesive soils, and test the effects of different tensile st rength. As a limiting 
case o f a strengthless material we use the HSL for water (IHolsapple and Housenl . 



20071 ). These equations read: 



gd \ "2+77 / p 



2v 
2+m 



respectively for cohesive soils and water. D is the crater diameter, v± is 
the normal component of the impact velocity, 5 is the impactor density, k, //, v 
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depend on the material and are derived from experiments. Their numerical 
values are k = 1.03, 1.17, fi = 0.41,0.55, respective l y for Eg. [2] and Eq. El 
while v = 0.4 in all cases ( Holsapple and Housenl . 2007 ). Concerning the 



strength, we may regard typical lunar regolith (Y ~ 10 5 dyne/cm 2 ) and bulk 
silicates (Y ~ 10 8 dyne/cm 2 ) as limiting cases for a silicate body. Highly 
under-dense (porous), aggregate materials having Y < 10 5 dyne/cm 2 can 
be ruled out because of Steins' stony composition (E-type). Concerning the 
higher limit, a more realistic estimate of the strength for an asteroid may 
be obtained considering that the strength depends on the asteroid size R, 
with larger bodies being weaker than smalle r ones of similar com position. 
Assuming that the strength scales as R^ 1 ^ 3 ( Asphaug et all 1996 ). Steins' 



hard rock strength may be as low as a few Y ~ 10 6 dyne/cm 2 . In the 
light of previous reasoning, we restrict the following analysis to the HSL 
for cohesive soils (for two representative values of tensile strength, namely 
Y = 10 5 , 10 6 dyne/cm 2 ). In fig. |5]the HSL obtained for different parameters 
are reported. 



A different approach is that proposed by lNolan et al.l ( 119961 ). They estimated 
the cratering scaling law (NSL, hereinafter) using hydrocode simulations. 
Their main result is the discovery of the so-called fracture regime, which 
occurs in between the two extreme situations represented by the strength 
and gravity regimes. Basically, small craters are formed in the classical way, 
with their size being controlled by the local strength. In large craters, on 
the other hand, the shock wave propagates ahead of the excavation flow, 
and therefore the material is totally fractured prior to its removal. If the 
amount of excavated material is large enough, the size of the resulting crater 
is controlled by the gravity. A NSL has been derived for Gas pra, which can 



be res caled to other asteroids using the approach proposed by lO'Brien et al 



( 120061 ). Note that in this scenario bodies can survive much larger impacts 



than predicted by HSL. NSL can be arranged in the following manner: 

D = «,»(« )'(»)' (4) 

where c, a, (3 and 5 have different values for the strength, fracture and 
gravity regimes. The transitions between strength to fracture regime and 
from fracture to gravity regime occur respectively at: 

( v \ 7 

d 'f = d o{—) ( 5 ) 
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where do, ^o and go are parameters computed for Gaspra. Numerical co- 
efficient used here (expressed in c.g.s. units) are^l: c = 35, 26.61, 161.4; a = 
1, 1.159, 0.78; = 0.56, 0.65, 0.44; 5 = 0,0, -0.22 respectively for the strength, 
fracture and gravity regimes, while 7 = 0= —0.56; 9 = —0.58 for all cases. 
Finally, d Q = 560 cm, d\ = 2.56 • 10 4 cm, vq = 5.0 km/s go = 0.448 cm/s 2 . 
Figure [5] shows the NSL rescaled to Steins. 

By comparing the SLs reported in fig. |5l a large degree of variation emerges. 
For a fixed impactor size d, the resulting crater diameter D may vary by 
more than a factor of 10. However, if we restrict ourselves to the observed 
crater size range on Steins and to the most likely scaling laws (NSL and HSL 
for 10 5 < Y < 10 6 dyne/cm 2 ) the variation is within a factor of 3. The 
difference can be partly explained by the fact that NSL overesti mates crater 



ing e fficiency since it neglects the shear resistance of materials (INolan et al. 
1996h . 



5. The model production function 

Using the considerations described in previous sections, it is possible to 
compute the differential distribution (<&(£))) of the number of craters with 
respect to their diameters expressed per unit time and surface area. The 
MPF can be obtained by: 

MPF(D) = / $(D)dD (7) 
Jd 

The distribution of craters for a given age t is simply obtained by: MPF(D) x 
t. These equations implicitly assume that all craters accumulate over time 
without interfering with previously formed craters (i.e. no crater erasing) 



2 Notice that the parameters involved in equations I4I5I6I are not totally independent 
because of continuity conditions at the boundary between two adjacent regimes. They 
can be written in the following way: /3 — 2£/(l — £),2a/£/(l — £),2£;5 = 0, 0,— £ for 
the strength, fracture and gravity regimes. Moreover, 7 = — 2£/(l — £);</> = 2£(1 — £ — 
a/)/[(l — £)(a/ — ct g )}; 9 = — £(l/(a/ — o> g )), where the subscript f and g stands for the 
fracture and gra vity regime and £ is a parameter that depends on the material and has 
been set to 0.22 ( O'Brien et all l2006f ). Note that the above equation have been simplified 
using a s = 1. 
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and that the flux is constant over time. The latter assumpt ion, according 



to lu nar chronology, is valid for ages less than ~ 3.7 Ga (e.g. iMarchi et al. 
20091 ). In figure |6] (upper panel) the fit of the MPF to Steins crater counting 
data is shown. 

The first important result is that the shape of the cumulative distribution 
cannot be satisfyingly fitted in the available size range. In particular, when 
fitting craters larger than ~ 0.6 km, the smaller craters are strongly under- 
represented. Similar kinks in the cumulative distributi on have been observed 
elsewhere, in particula r on the Moon and Mars (e.g. iHiesinger et al.l . 12002 ; 
Haruyama et all 120091 ). and usually are attributed to episodes of crater eras- 
ing which are more effective for small diameters. A similar lack of small 
craters, but for sizes < 10 m, has also recently been observed on Itokawa 

In fig. [6] (upper panel) the best 



„ ~„ 

(Hirata et al.. 


2009; 


Michel et al.. 


2009) 



fit for large crater diameters is shown. The corresponding age, derived by 
using the NSL is 0.1 ± 0.02 Ga. This estimate is likely to represent a lower 
limit, since the NSL n eglects shear resist ance and therefore tends to overesti- 
mate the crater sizes (INolan et al.l . Il996l ). The HSL derived age ranges from 
0.28 ± 0.06 Ga to 0.94 ± 0.2 Ga, for Y = 10 5 , 10 6 dyne/cm 2 , respectively. 
Model ages are determined through a x 2 fitting. Formal errors are estimated 
considering a variation of ±30% around the minimum \ 2 . The cratering pro- 
cess may be more complicated than assumed so far, in particular the MPF 
may vary over time. The main reason for such time dependence is that craters 
may erase over time. A number of processes responsible for crater erasing 
on small bodies have been indentified: local and global jolting , cumulative 
seismic shaking and su perposition of craters ( Greenberg et al. . 19941 1996; 



Richardson et al 



2006) 



20041 ) . Such effects can be modeled (I O'Brien et al 
and t he resulting MPF can be written in the following manner (IMarchi et al 
2009( 1: 



MPF (D,t) 



<$>(D,t)£{D,t) dD 



(8) 



D 



where the function £(D,t) is the ratio of the final number of craters, 
erasing included, to the total number (i.e. erasing excluded). The mentioned 
erasing process depends on several parameters. As f or the regolith j olting 
and superposition, we use the parameters adopted in lO'Brien et al.l (120061 ) 
for Gaspra. Cumulative seismic shaking has not been applied in our simula- 
tions due to the lack of detailed information on regolith mobility on Steins. 
Figure El (lower panel) shows the effects of crater erasing on Steins' age de- 
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termination. Also in this case, even though the MPF is shallower, the small 
craters are not accurately fitted by the models. Derived ages are now in- 
creased with respect to what was found neglecting crater erasing. Focusing 
again on D > 0.6 km, we derive an age of 0.154 ± 0.035 Ga for the NSL. In 
the case of the HSL, the age becomes 0.49 ± 0.18 Ga and 1.6 ± 0.5 Ga for 
Y = 10 5 , 10 6 dyne/cm 2 , respectively. 

We also investigated whether or not the observed crater population is satu- 
rated, i.e. has reached an equilibrium point where new craters erase old ones 
leaving unchanged the overall distribution. We find that the crater popu- 
lation has not yet reached saturation. To see this point, in figure |6] (lower 
panel) we overplot the MPF for 0.5 Ga (NSL) and 3.6 Ga (HSL). These 
ages have been chosen in order to have the corresponding MPFs above all 
the cumulative data points. Both MPFs are clearly separated by the best-fit 
curves, indicating that the saturation is not reached yet, therefore age as- 
sessment is possible. 

It is interesting that also when taking into account the erasing of craters 
by superposition and regolith jolt, smaller craters are still strongly under- 
represented. A number of possible explanations for the origin of this kink 
may be invoked. For instance, small craters may have been erased by re- 
golith displacement due to the effect of cumulative seismic shaking. This 
process has been demons trated to expla in the deficit of sma ll craters on Eros 
( jRichardson et al.l . I2004J ) and Itokawa (IMichel et al.l . 120091 ) . In order to ad- 
dress this issue, we show the Steins' crater distribution on a R-plot (see figure 
[7]). R-plots are useful in order to analyze fine details of crater distributions. 
Data from Itokawa and Eros are also overplotted. The comparison among 
Itokawa, Eros and Steins is interesting mainly for the purpose of the small 
craters' depletion, although, it is not so straightforward because of the dif- 
ferent size range of craters. However, some interesting comments can still be 
drawn. 

Itokawa and Eros show a marked depletion of small crater, for D < 0.1 km. 
This fact has been interpreted as the result of crater erasi ng triggered by cu- 



mula t ive seismic shaking due to repeated small impacts (jRichardson et al. 



2004 ; IMichel et al.l . 120091 ) . This process is able to reproduce accurately the 



observed trend on Itokawa and Eros (indicated for clarity by the line l in 

fig. CD- 

As for Steins, the crater distribution for D < 0.5 — 0.6 km shows an overall 
similar behavior, indicating that also on Steins the cumulative seismic shak- 
ing may be a viable explanation for the lack of small craters. Nevertheless 
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there are some fine details that differ from Itokawa and Eros which are worth 
to be analyzed. To better show these discrepancies, in fig. [7] we draw the line 
li which is parallel to 1$. Despite the large error bars, the 3 left-most bins 
(D < 0.3 — 0.4 km) exhibit a steeper trend than that of Itokawa and Eros, 
although they might be compatible due to the relatively large statistical er- 
rors. This may be an indication that the erasing on Steins had a different 
origin, or at least that the cumulative shaking is not the unique responsible 
of the observed lack of small craters. We suggest that the Steins observed 
peculiar crater distribution could retain the footprint of an intense episode of 
erasing triggered be a single event, likely the formation of Ruby, which would 
have erased preferentially craters below D ~ 0.3 — 0.4 km. Note that the size 
bin corresponding to the Ruby crater in the cumulative distribution is close 
to the best fit model curve (see fig. [6]). This means that the formation of 
such large craters is already accounted for by the model and yet this is not 
sufficient for explaining the lack of small craters, at least with the parameters 
used in this study. Therefore, we may argue that the formation of the Ruby 
crater happened in a more recent time than the best fit age. 
An alternative explanation is connected to the YORP evolution that may 
have triggered regolit h mobility and efficiently erased small craters (e.g. see 
Scheeres et al. , 2007t ). 



Without detailed knowledge of the internal nature of Steins and its regolith 
thickness it is very difficult to draw a firm conclusion. Nevertheless, erasing 
triggered by the Ruby impact seems more likely than the one induced by 
the YORP effect. This conclusion seems at least partially confirmed by the 
fresh appearance of the Ruby crater, with its sharp rims and high depth-to- 
diameter ratio. 

In any case, we may use the crater distribution to constrain the epoch of 
the putative "impulsive erasing". In order to correctly achieve this result 
a detailed knowledge of the erasing process is required. This is however far 
beyond the scopes of the present paper, hence we limit our discussion consid- 
ering two cases. First of all, the case where all small craters (D < 0.35 km) 
were erased at the same time. Through MPF fitting of the observed cumu- 
lative distribution for small craters, the approximate time lapsed since the 
erasing can be derived. Considering only the craters for D < 0.35 km, for 
NSL we obtain an age of 0.032 ± 0.004 Ga, which becomes 0.072 ± 0.01 Gy 
and 0.237 ± 0.03 Gy, for HSL and Y = 10 5 , 10 6 dyne/cm 2 , respectively (see 
fig. O lower panel). Note, however, that the accumulations of craters during 
this time would have produced a steeper slope than the observed one result- 
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ing from face value of data points. Nevertheless, the statistical errors are too 
large to conclude whether or not the model distribution is really different 
to the observed trend. In this respect, it cannot be ruled out that seismic 
shaking might be entirely responsible for the slope of the crater size distri- 
bution in that size range, in this case the observed trend cannot be used to 
constrain the timescale of Ruby's formation. 

A lower limit for the epoch of the impulsive erasing can be derived consid- 
ering that it erased most or all of the smallest craters (D ~ 0.2 — 0.3 km), 
but only some of the larger craters in order to reproduce the observed slope. 
In this case, the timescale could be set by the timescale to form just the 
smallest craters. The time required to accumulate the 4 observed craters 
having D = 0.2 - 0.3 km spans from ~ 2 Ma to ~ 10 Ma, for NSL and HSLs 
respectively. 

It is possible, using Poisson statistics, to compute the expected probabilities 
that the impulsive erasing due to Ruby formation occurred at the computed 
times. Assuming that during the crater retention ages the formation of only a 
single crater of the size of Ruby occurred, the probabilities that Ruby event 
occurred in the estimated recent times are: from ~1 to ~17% (NSL) and 
from ~0.6 to ~12% (HSL), respectively for the two different estimates of the 
impulsive event times reported above. 

Another consequence of the formation of the Ruby crater would be the mix- 
ing of the regolith layer, with subsequent reset of the optical properties. It 
is interesting that detailed investigatio ns acro s s the surface of Steins have 



shown little or no spectral variability (ILeyratl . 120101 ) . This is somehow in 



contrast with other asteroids visited by spac ecraft, all showin g a certain de- 
gree of alteration due to space weathering ( IChapmanl . 120041 ) . The lack of 
spectral variability across Steins can be explained in different ways, although 
it must be noticed that space weathering response on E-type materials has 
not been investigated in details yet. A possibilty is that the lack of spec- 
tral variability may be an indication of Steins' insensitive response to space 
weathering alteration due to its specific composition. However, as demon- 
strated by lLazzarin et al.l ( 120061 ). all asteroidal spectral types are found to 
show some degree of alteration. In this case, Steins would have either a fresh 
or a totally saturated surface. The fresh surface scenario seems more likely 
given the fact that there is no spectral variation in proximity of small -and 
therefore relativey young- craters. This conclusion would be in agreement 
with the relatively young age estimated for the formation of the Ruby crater. 
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6. Conclusions 



In this work we used crater counting, morphological analysis and im- 
pactor population modeling to constrain Steins' cratering retention age. The 
derived ages vary according to the SLs used. In particular, using NSL and 
crater erasing the derived age is 0.154±0.035 Ga, while using HSL and crater 
erasing the age ranges from 0.49 to 1.6 Ga, according to the values of strength 
used. Moreover, the modeling of the crater erasing processes shows that the 
observed crater density is not saturated (at least for the parameters adopted 
here). The mean co llisional age of Steins is estimated to be ~ 2.2 Ga (e.g. 



Marchi et al.l . 120061 ) . Interestingly, similar numbers apply also for Gaspra. 



Analogue conclusions might be also val id for the near-Ear th objects Eros, 
whose mean collisional age is ~ 1.7 Ga (IMarchi et al.l . 120061 ) while cratering 



age using NS L give s 0.12 Ga flO'Brien et all . 120061 ) or, using HSL, 1-2 Ga 



( jMichel et al.l . 120091 ) . The larger bodies Ida and Mathilde have crater pop- 
ulations either close to the saturation or saturated, and consequently their 
cratering age estimate is less constrained. Despite the low number statistics, 
the cratering age of the main belt asteroids smaller than ~ 20 km seems to be 
systematically younger than their collisional age. This result, if confirmed by 
further studies, would have important implications on main belt collisional 
models. 

Notably, the shape of Steins crater size distribution shows a kink for diame- 
ters smaller than 0.5-0.6 km, which may require a recent episode of intense 
erasing, although seismic shaking could have potentially played a role in pro- 
ducing the observed distribution as well. We also attempt to constrain the 
epoch of such episode, possibly associated to the Ruby crater formation. Fo- 
cusing on the small diameter end (D < 0.35 km) of the crater cumulative dis- 
tribution, and adopting the MPF fitting we obtain an age of 0.032 ±0.004 Ga 
using NSL, from 0.072 to 0.237 Ga using HSL. Under the assumption that 
the formation of the Ruby crater erased all craters < 0.5 — 0.6 km, the above 
ages could possibly indicate the time since the occurrence of the Ruby event. 
A lower limit to the age of this event, can be derived considering the time 
required to accumulate the observed craters in the range 0.2 — 0.3 km. This 
produces an age from ~ 2 Ma to ~ 10 Ma, for NSL and HSLs respectively. 
The derived ages vary up to a factor of ten depending on the SL and the 
tensile strength used. In the present work we investigated the effects of a 
relatively large range of SLs and Y. However, in the light of our global un- 
derstanding of Stein properties, we favor a crater retention age ranging from 
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~ 0.15 to ~ 0.5 Ga, and a kink related event that could be as young as a 
few Ma up to some tens of Ma. 

As a final remark, note that the conclusions derived in this paper are based 
on the bona fide crater distribution. The general scenario outlined (age, 
depletion of small craters) remains valid also if considering all crater-like fea- 
tures (see fig. El lower panel). In the latter case, however, the age of the 
reset event is about a factor of two higher. 
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Figure 1: WAC image (left) at the close approach. The image was acquired with a narrow 
passband filter with central wavelength of 295.9 nm and FWHM=10.9 nm. The spatial 
resolution is 80 m/px. Since craters close to the limb are not easily detectable, the actual 
region used for crater count is reduced by one pixel around the limb. The pit-chain is 
clearly visible close to the right-hand terminator. The NAC image (right) shows the 
large depression which is approximately in the antipodal position of the alignment of pits 
and Ruby crater. The image was acquired with a narrow passband filter with central 
wavelength of 805.3 nm and FWHM=84.5 nm. The spatial resolution is 98 m/px. 



19 



0.1 



0.01 



All crater-like features 
Bona fide craters 



_L 



0.1 



Crater diameter (km) 



Figure 2: Cumulative distribution of all crater-like features and bona fide craters. The 
latter is the distribution used in the res t of th e paper for age assessment. Both distributions 
contain features > 3 pixels (see iBesse , 20ld for more details). Notice that the bona fide 
craters distribution does not contain the depressions forming the chain-like features (see 
text for details) and a couple of uncertain craters close to the Ruby's rim. The resulting 
distribution shows a remarkable paucity of small (D < 0.5—0.6 km) craters when compared 
to the distribution of all detected crater-like features. Error bars are estimated on the basis 
of Poisson statistics of counts. 
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Figure 3: Steins' cumulative impactor size distribution (upper panel) and impact velocity 
distribution (lower panel) us ed in the present wo rk. For a comparison, the lunar impact 
distributions are also shown ( Marchi et al. . 20091) . 
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Figure 4: The ratio of the diameter of the biggest crater divided by the average diameter 
of various asteroids (triangles, left ordinate) and the third root of the specific energy 
for disruption (dia monds) and shattering (ast erisks) of a basaltic body impacted with a 
velocity of 5 km/s ( Benz and AsphaugL 19991 right ordinate). The relation between the 
left and right ordinate is arbitrary. The critical impact is expected to be between the 
shattering and the disruption limit. For Steins we used an average diameter of 5.3 km, 
a size of the largest crater of 2 km and assumed a density (nee ded for the ener gy of 
disruption) of 2 g/cni 3 . Data for the other asteroids are taken from lAsphau g (2008). For 
Amalthea and Vesta the specific energy of disruption is above the plot range. The figure 
shows that the ratio between size of the largest crater and body size is not particularly 
high on Steins. However, since the energy needed to shatter a body is at a minimum 
at approximately the size of Steins, an impactor large enough to create the Ruby crater 
wo uld be expe c ted to shatter (but not necessarily disrupt) Steins. The model calculations 
bv I Jutzi et al. ( 2010l ) suggest that this actually may have been the case: In many cases 
an impactor of the size needed to create the Ruby crater will strongly damage the original 
body, but Steins' global shape will remain approximately the same and after the event the 
crater will be clearly visible. 
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Figure 5: Relationship between impactor diameter (d) and crater diameter (D) according 
to the scaling laws discussed in this work. We also indicated the approximate impactor 
diameter for catastrophic disrupt (e? C£ j) of Steins, obtained for an unfractured silicate body 
and HSL. The limiting case of the water scaling law (not used for dating purposes) is 
shown only for a comparison. The gray area indicates the crater size range detected on 
Steins. 
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Figure 6: Upper panel: Steins' age estimates obtained with the MPF and no crater erasing. 
The best fits have been performed considering only craters for D > 0.6 km. Notice the 
strong lack of observed craters for D < 0.6 km in comparison to the models. Lower panel: 
Steins' age estimates obtained with the MPF and crater erasing. The best fits for small 
diameters (D < 0.35 km) is also shown, possibly indicating the time of the formation of 
the Ruby crater (see text for further details). For a comparison, the distribution of all 
crater- like features is also shown (errors bars are not shown for simplicity, see Fig. [2j. 
To adreess the crater saturation issue a couple of MPFs for older ages have been also 
overplotted (see text). 
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Figure 7: Steins' crater data shown on a Relative plot (R-p l ot). R- values are com- 
puted according to Crater Analysis Techniques Working Group ( 1979h . For comparison, 
Itokawa and Eros d ata are also shown. Itokaw a data hav e been computed from table 1 of 
Hirata et al. (2009). while Eros data are from IChapmanl (|2002h . Dashed line l Q indicates 



the average slopes for the left- most part of the Itokawa/Eros distributions (D < 0.15 km), 
showing a paucity of small craters. For Itokawa an d Eros this paucity of small craters has 
been interpreted as the result of seismic shaking (jRichardson et al. . l2004t iMichef et all 
2009). Steins also shows a similar behavior (see line l\ which is parallel to l a ) but for 



D < 0.35 km the trend seems to have a steeper slope than Itokawa/Eros, possibly due to 
different processes of crater erasing (see text for details). 
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